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ABSTRACT 

We propose to study cosmic reionization using absorption line spectra of high-redshift Gamma 
Ray Burst (GRB) afterglows. We show that the statistics of the dark portions (gaps) in 
ORB absorption spectra represent exquisite tools to discriminate among different reioniza- 
tion models. We then compute the probability to find the largest gap in a given width range 
[Wmax, Winax + dW] at a flux threshold Fth for burst afterglows at redshifts 6.3 < z < 6.7. 
We show that different reionization scenarios populate the (M^maxj Fth) plane in a very dif- 
ferent way, allowing to distinguish among different reionization histories. We provide here 
useful plots that allow a very simple and direct comparison between observations and model 
results. Finally, we apply our methods to GRB 050904 detected at z = 6.29. We show 
that the observation of this burst strongly favors reionization models which predict a highly 
ionized intergalactic medium at z 6, with an estimated mean neutral hydrogen fraction 
xm = 6.4 ± 0.3 X IQ-^ along the hne of sight towards GRB 050904. 

Key words: cosmology: large-scale structure of Universe - intergalactic medium - gamma- 
ray: bursts 
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1 INTRODUCTION 

In the last few years, our knowledge of the high- 2: Universe and 
in particular of the reionization process has been enormously in- 
creased mainly owing to the observation of quasars by the SDSS 
survey (Fan 2006) and CMB data (Hinshaw et al. 2007, Page et 
al 2007). Long gamma ray bursts (GRB) may constitute a comple- 
mentary way to study the reionization process avoiding the proxim- 
ity effects and possibly probing even larger redshifts. This has now 
become clear after the detection of five GRBs at 2 ^ 5, over a sam- 
ple of about 200 objects observed with the Swift satellite (Gehrels 
et al. 2004). The current record holder is GRB 050904 at z = 6.29 
(Tagliaferri et al. 2005, Kawai et al. 2006). Totani et al. (2006) have 
used this object to constrain the ionization state of the intergalactic 
medium (IGM) at high redshift by modeling its optical afterglow 
spectrum. They report the evidence that the IGM was largely ion- 
ized already at z — 6.3. The best-fit neutral hydrogen fraction is 
consistent with zero with upper limit xm < 0.17 (< 0.6) at 68% 
(95%) C.L. 

Several authors (Natarajan et al. 2005, Daigne, Rossi & 
Mochkovitch 2006, Bromm & Loeb 2006, Salvaterra et al. 2007a) 
have computed the number of high-2 GRBs detectable by Swift. In 
spite of model details, all these different studies consistently predict 
that a non-negligible fraction (up to ~ 10%) of all observed GRBs 
should lie at very high redshift. On the basis of these results, several 
GRBs at 2 > 6 should be observed by Swift in the near future. De- 



pending on the Swift/BAI trigger sensitivity and on model details, 
Salvaterra & Chincarini (2007) found that ~ 2 - 8 GRBs can be 
detected above this redshift during every year of mission. From the 
observational point of view, Ruiz-Velasco et al. (2008) find that the 
fraction of Swift GRBs at 2 > 6 should not exceed the conserva- 
tive upper limit of 19%. Moreover, future X- and Gamma-ray mis- 
sions will increase rapidly the sample of high-2 GRBs, possibly up 
to 2 ~ 10 (Salvaterra et al. 2007b). Once detected, spectroscopic 
follow-up observations of high-2 GRBs require a rapid trigger of 8- 
meter, ground based telescopes. This can be done by pre-selecting 
reliable candidates on the basis of some promptly- available infor- 
mation provided by Swift, such as burst duration, photon flux, the 
lack of detection in the UVOT l^-band, and the low Galactic ex- 
tinction (Campana et al. 2007, Salvaterra et al. 2007a). We note that 
cosmological time dilation helps keeping the flux bright, since ob- 
servations will sample the afterglow early phases, even a few days 
after trigger. 

The spectra of high-redshift sources (as QSOs and GRBs) 
bluewards of the Lya are characterized by dark portions (gaps) pro- 
duced by intervening neutral hydrogen along the line of sight. The 
use of various gap statistics in QSO spectra has been recently rec- 
ognized as a very powerful tool to constrain the IGM ionization 
state (Paschos & Norman 2005; Fan et al. 2006; Gallerani, Choud- 
hury & Ferrara 2006, hereafter G06; Gallerani et al. 2007, hereafter 
G07). For example, by comparing the statistics of these spectral 
features in a sample of 17 observed QSOs with Lya forest simu- 
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lations, G07 concluded that the HI fraction, xm, evolves smoothly 
from 10"*'* at z = 5.3 to 10~^'^ at z — 5.6, with a robust up- 
per limit xm < 0.36 at z — 6.3. These results encourage the 
application of such analysis to GRBs. There are several advantages 
promised by such an attempt. First, GRBs are soon expected to be 
found at redshifts higher than those typical of QSOs; second, and 
contrary to the massive hosts of QSOs, they reside in "average" cos- 
mic regions only marginally affected by local ionization effects or 
strong clustering; finally, they are bright and their afterglow spec- 
tra closely follows a power-law, making continuum determination 
much easier 

In agreement with the WMAP3 results (Spergel et al. 2007) 
we assume a flat universe with Q,m ~ 0.24, Qa ~ 0.76, f^^/i^ = 
0.022, h — 0.73. The parameters defining the linear dark matter 
power spectrum are n — 0.95, dn/dlnfc = 0, as = 0.82. Mpc 
are physical unless differently stated. 



2 GRB EMISSION PROPERTIES 

We have built a synthetic GRB afterglow emission spectrum start- 
ing from the observed spectral energy distribution and time evolu- 
tion of the most distant GRB detected up-to-now, i.e. GRB 050904 
at z = 6.29 (Tagliaferri et al. 2005, Kawai et al. 2006). The unab- 
sorbed afterglow spectrum of GRB 050904 can be parametrized as 
F{iy) oc u"t'^, with a = ~1.25to[ll (Haislip et al. 2006, Taglia- 
ferri et al. 2005). The observed temporal decay up to 0.5 days from 
burst is well described by a power-law index of —1.36 ± 0.06, fol- 
lowed by a plateau phase with f3 — —0.82 ± 0.15 lasting until 2.6 
days after burst (Haislip et al. 2006). The further afterglow evolu- 
tion can be described assuming (3 = —2.4 ± 0.4 (Tagliaferri et 
al. 2005). Finally, we normalize the intrinsic GRB 050904 optical 
spectrum in order to reproduce the flux of 18 fj,}y as measured 
at 1 day from burst in the J band (Haislip et al. 2006). Although we 
have experimented with other choices to explore the model sensi- 
tivity to redshift, we will mostly show results for GRBs located in 
the range 6.3 < z < 6.7, which appears to be crucial to properly 
follow the latest (i.e. overlap) phases of reionization history. From 
the observed afterglow spectral evolution of GRB 050904, we com- 
pute the rest frame spectrum between Lya (1215.67 A) and Ly/3 
(1025.72 A). 



3 ABSORPTION SPECTRA 

The ultraviolet radiation emitted by a GRB can suffer resonant Lya 
scattering as it propagates through the intergalactic neutral hydro- 
gen. In this process, photons are removed from the line of sight 
(LOS) resulting in an attenuation of the source flux, the so-called 
Gunn-Peterson (GP) effect. Thus, the observed flux Fobs is given 
by -Fobs ~ F(p) e~'"'^P, where top is the GP optical depth. To 
simulate the tgp distribution we use the method described by G06 
and further revised in G07, whose main features are summarized 
as follows. Mildly non-linear density fluctuations giving raise to 
spectral absorption features in the intergalactic medium (IGM) are 
described by a Log-Normal distribution, which has been shown to 
fit observational data at redshifts 1.7 < z < 5.8 (Becker et al. 
2007). For a given IGM equation of state (i.e. temperature-density 
relation), the mean HI fraction, xhi, can be computed from pho- 
toionization equilibrium (an excellent approximation under the pre- 
vailing physical conditions) as a function of the baryonic over- 
density, A = p/p, and photoionization rate, F, due to the ultra- 
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Figure 1. Synthetic absorption spectra at different times after the burst 
(t=0.1, 1.0, 3.4 days from the bottom to the top) of a GRB at 2 = 6.3, for 
the ERM. The thin black horizontal lines denote the flux threshold F°^^ 
used to define gaps. The thick dashed cyan lines show the largest dark gap 
along each lines of sight. 

violet background radiation field. These quantities must be deter- 
mined from a combination of theory and observations; here we fol- 
low the approach of Choudhury & Ferrara (2006, CF06 hereafter). 
Their study allows for three types of ultraviolet ionizing photons: 
QSOs, PopII and PopIII stars. As the contribution of PopIII stars 
is negligible for the redshift of interest here, we neglect this type 
of sources. The CF06 model contains two free parameters: (i) the 
star-formation efficiency /* , and (ii) the escape fraction /esc of ion- 
izing photons from galaxies. These are calibrated by a maximum- 
likelihood procedure to a broad data set, including the redshift evo- 
lution of Lyman-limit systems, Lya and Ly/3 GP optical depths, 
electron scattering optical depth, cosmic star formation history, and 
number counts of high redshift sources. 

Currently, the available data can be explained by two dif- 
ferent reionization histories, corresponding to different choices 
of the free parameters: (i) an Early Reionization Model (ERM) 
(/* = 0.1;/esc = 0.07), and (ii) a Late Reionization Model (LRM) 
(/* = 0.08; /esc = 0.04). These two models bracket the range of 
possible reionization histories. For example models in which reion- 
ization occurs earlier (later) than in the ERM (LRM) would over- 
produce the electron scattering (GP) optical depth. The global prop- 
erties of the two reionization models are shown in Fig. 1 of G07 so 
we do not report them here. It suffices to say that in the ERM the 
volume filling factor of ionized regions, Qhii = Veii/Kot = 1 
at z < 7\ in the LRM Qhii evolves from 0.65 to unity in the red- 
shift range 7.0-6.0, implying that the reionization is still in the pre- 
overlap stage at z > 6. Both ERM and LRM by construction pro- 
vide an excellent fit to the mean neutral hydrogen fraction evolution 
experimentally deduced from the GP test. Having selected all the 
necessary physical evolutionary properties of the most likely reion- 
ization histories, we can now describe how we plan to use GRBs to 
discriminate between them. 



4 REIONIZATION TESTS 

The main idea we propose in this paper is to exploit the statistics of 
the transmissivity gaps imprinted by the intervening IGM neutral 
hydrogen on the otherwise smooth power-law spectrum of high- 
redshift GRBs. On general grounds, we expect that at any given 
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redshift, but particularly above z = 6, where differences become 
more marked, the value of xhi is higher in the LRM than in the 
ERM. As a result, if reionization completes later, absorption spec- 
tra are characterized by wider and more numerous gaps, defined as 
contiguous regions of the spectrum with observed flux lower than 
a given flux threshold over a rest frame wavelength interval larger 
than 1 A. Moreover, the level of neutral hydrogen fluctuations along 
the LOS to the GRB differs in the two models: in principle, the 
GRB flux decay could be used as a tunable high-pass filter which 
allows to study the growth of gaps with time. As the time after 
the burst increases, the gaps become larger. In fact, the progres- 
sive fading of the unabsorbed afterglow produces a corresponding 
attenuation of the observed flux, as can be seen from Fig. 1. As 
the flux falls below the flux threshold F°^^ (thin black horizon- 
tal lines) the transmission windows disappear, causing the gaps to 
expand. The thick dashed cyan line shows the largest gap (M^max) 
along the LOS; the largest gap width (LGW) increases as the GRB 
flux fades. However, in practice, the analysis of the LGW time evo- 
lution can be equivalently applied at a fixed time after the burst to 
the transmitted flux Ftransm = Fohs/F(i/) varying the threshold^ 

We have derived the evolution of the LGW found in S3mthetic 
absorption spectra varying the flux threshold used to define gaps. 
Suppose that a GRB is observed at a given redshift zgrb . We can 
then ask what is the probability that in its afterglow spectrum the 
largest gap, defined by Fth, is found within a given width range^ 
[M^max , Wmax . We cau draw this probability from our model 
by analyzing a large number of synthetic spectra corresponding to 
different GRB redshifts «grb and counting the number of LOS 
having the largest gap in the given range, varying the Fth between 
and L Clearly, in real cases, a minimum threshold Fth,min is set 
by the noise associated to the spectrum. 

We have computed this statistics for afterglow spectra of 
GRBs with 6.3 < «grb < 6.7. In fact, according to the most 
recent estimates involving the lower trigger sensitivity, the Swift 
mission is expected to increase the sample of high-^; GRBs, possi- 
bly up to 2: ~ 10 (Salvaterra et al. 2007b). 

We show the resulting largest gap probability isocontours in 
Fig. 2 for ZQRB = 6.3 (top panels), «grb = 6.5 (middle panels), 
and ^GRB = 6.7 (bottom panels); the left (right) panels refer to 
the ERM (LRM) case. We use 300 LOS in order to ensure statis- 
tical convergence of the results. The isocontours correspond to a 
probability of 5%, 10%, and 40%. The differences caused by the 
two different reionization histories are striking, since the ERM and 
LRM populate the (Wmax, fth) plane in a very different way, as 
shown in Fig. 2. For example, at zcjrb = 6.3, the probability to 
find the largest gap with M^max = 50 A for Fth = 0.1, is « 40%, 
in the ERM, while it is only of the « 5% in the LRM, being most 
of the largest gaps wider than 50 A. 

Moreover, in the LRM the width of the gaps is a factor « 2 
times wider than in the ERM. In particular, for zqrb = 6.5 and 
Fth = 0.2, we find that in the ERM 40 < W^max < 170 A, while 
in the LRM gaps are mostly larger than 80 A, for the same value of 
the flux threshold. 

Finally, in the ERM the LOS becomes completely dark for flux 
threshold values higher than in the LRM case, as a consequence of 
the higher transmitted flux predicted by the former. For example, at 



Note that we define F^^'^ as the threshold which defines gaps in the ob- 
served flux and Fth in reference to the transmitted one. 
^ In this work we assume dW = 20A. 



2 = 6.7, in the ERM (LRM) the gap size is equal to the simulated 
LOS length (« 190A) for Fth « 0.4 (0.15). 

The results point also towards two very important, and poten- 
tially interesting, redshift-dependent features, both related to the 
thickening of the hya forest, which becomes on average more 
neutral towards higher redshifts. First, there is an overall shift to- 
wards larger Wmax values of both ERM and LRM curves from 
^GRB = 6.3 to ^GRB = 6.7. Second, the rate at which the gap 
width grows and saturates with Fth strongly correlates with red- 
shift. 

Thus, the analysis of the largest gap evolution with the flux 
threshold shows that it is is easier to discriminate among different 
reionization histories using the highest redshift available sample of 

GRBs. 

Fig. 2 allows a straightforward comparison between data and 
model results. In practice, iirformation on the reionization history 
from a « > 6 GRB can be easily extracted through the following 
the steps: (i) from the J-band afterglow flux, Fj, extrapolate the 
continuum as F(i^) — Fj {v / vj)" ; (ii) compute the transmitted 
flux Ft ransm = Fobs/F(i/); (iii) Determine the LGW in Ft ransm, 
varyingFthin [Fth.min; 1], (iv) compare the result with Fig. 2. This 
simple procedure will allow to dicriminate between different reion- 
ization scenarios. 



5 THE CASE OF GRB 050904 

We apply the LGW analysis to the observed flux of the GRB 
050904, whose spectrum has been obtained 3.4 days after the burst 
(Kawai et al. 2006). For a fair comparison with observations, simu- 
lated spectra have been convolved with a gaussian of FWHM = 300 
km/sec, providing R = 1000, and then rebinned to a resolution of 
R = 600. Finally, we add noise to the simulated data such that the 
flux F in each pixel is replaced by F + G(l) an, where G(l) is a 
Gaussian random deviate with zero mean and unit variance, and 
an ~ 0.05 is almost the mean noise r.m.s deviation, observed in 
the GRB 050904 spectrum by Kawai et al. 2006. 

In this work the lower flux threshold value that we use is 
-Flt'min = 0.3 X 10-^*erg s-^cm-^A-^ which is > 3-ct„. Such 
conservative choice of Fth''" excludes the possibility that the widths 
of the gaps are underestimated because of some spurious transmis- 
sion peaks. For Ft°h' = (0.30; 0.35; 0.40) x 10-^*erg/s/cmVA, 
the observed largest dark gaps sizes are Wmax = (51; 78; 141) A, 
respectively, in the GRB rest frame. These values refer to the dark 
regions extending in the observed wavelength region ~ [7850 — 
8220]; [7850 - 8420]; [7850 - 8874] A, respectively. It is worth 
noting that the third wavelength range contains the spectral region 
immediately blueward the Lya emission, which could be affected 
by the presence of a DLA, as claimed by Kawai et al. (2006). Thus, 
the corresponding quoted size for the dark gap should be consider 
as an upper limit. 

The black filled circles in Fig. 3, corresponding to the val- 
ues of Wmax and F^^^ for GRB 050904, show that the ERM is 
favored by the data. In particular, the probability to find a largest 
gap of Wmax ~ 51 A, for Ft^' = 0.30 x 10"^*crg/s/cmVA is 
~ 40% if reionization occurred early, i.e. almost half of the LOS 
contains a largest gap of this size for a burst with the redshift and 
flux of GRB 050904. Such probability drops for a late reioniza- 
tion model to ^ 10% clearly indicating that in this scenario the 
GRB 050904 observation represents a much rarer event. Moreover, 
in the ERM, the probability to find a largest gap of Wmax w 78 A, 
for Ft°h " = 0.35 X 10""erg/s/cmVA is almost 3 times the LRM 
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Figure 2. Isocontours of the probability that the afterglow spectrum of a GRB at redshift zgrb, contains a largest gap of size in the range [Wmax, Wmax + 
dW], with dW = 20 Afor a flux threshold . The left (right) panel shows the results for the ERM (LRM). The top (middle, bottom) panel shows the results 
for 2GRB = 6.3 (6.5, 6.7). The isocontours correspond to probability of 5%, 10%, and 40%. 



one (27% vs 10%). Finally, for the third observed point the LRM 
appears to be more probable than the ERM by a factor 27:5. This 
value should be considered an upper limit, since the third point has 
been obtained analyzing a spectral region which could be affected 
by a DLA contamination (Kawai et al. 2006). For this reason, the 
corresponding gap size should be smaller than 14lA, thus being 
more consistent with ERM predictions, as the other observational 
data. Without taking into account the latter observed point, we find 
that the ERM is 1 1 times more probable than the LRM; if such data 
point is included, the above vantage factor reduces to 2. We have 
checked the impact of uncertainties on the continuum flux by vary- 
ing its slope within the quoted error on a (see Sec. 2). We have 
found that the results shown in Fig. 3 are negligibly affected. 
Our results are due to the fact that many, relatively large, opaque 
stretches are present in the IGM at z — 6.3 for the LRM, result- 
ing in an average large Wmax- In the case of ERM, reionization is 
almost complete already at 2 ~ 7 and smaller largest gaps are ex- 
pected at the redshift of GRB 050904. 

Although a large sample of high redshift GRBs is required before 
we conclude that a model in which reionization ended by 2 ~ 7 
is favored by the data, the discriminating power of the proposed 
method is already apparent. 

5.1 Neutral hydrogen measurements 

We have compared the gap statistics as predicted by two differ- 
ent reionization models. At z=6, the ERM predicts a mean neutral 
hydrogen fraction xm ~ 10~*, while in the LRM xm ~ 0.03. 
As shown in the previous Sections, the gap statistics is sensitive 



to the neutral hydrogen amount along the LOS to the background 
source. Thus, we can exploit the better agreement of the ERM pre- 
dictions with observed data to provide a measurement of the neu- 
tral hydrogen fraction xm along the LOS to the GRB 050904. To 
this aim, from the ERM sample used in the analysis, we select 
those LOS whose largest gap is consistent with observations. Thus, 
we derive the mean neutral hydrogen fraction xm along the syn- 
thetic LOS characterized by W^ax e [51 A ± dW] for Ft^= = 
0.30 X 10"'* erg s"' cm-^ A"' and W^a^ e [78 A ± dW] for 
F°^' = 0.35 X 10"'** erg s"' cm"^ A"'. We find that the ob- 
served LGW in the GRB 050904 afterglow spectrum are consistent 
with Xm = 6.4 ± 0.3 x lO"''. 



6 CONCLUSIONS 

We have proposed to investigate cosmic reionization using absorp- 
tion line spectra of high-z GRB afterglows. The evolution of the 
largest gap width (LGW) as a function of the flux threshold Fth 
used to define gaps shows marked differences in the two reion- 
ization models favored by present data. In particular, (i) the LGW 
is typically 2 times wider in the LRM than in the ERM, for a 
fixed Fth', (ii) the difference between ERM and LRM in terms of 
(W^max,-Fth) increases with z; (iii) an overall shift towards larger 
WVax values for a fixed Fth is found in both models towards higher 
redshifts. This analysis has shown that we can robustly distinguish 
among different reionization histories: improved results can be ob- 
tained if data is collected promptly after burst detection and by us- 
ing GRBs at the highest available redshifts. 
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Figure 3. Isocontours of the probability tliat the afterglow spectrum associated with a GRB at redshift ^grb = 6.3, contains a largest gap of size in the range 
[VFrnax, VVmax + dW], with dW = 20 A, for a flux threshold F^^^. The left (right) panel shows the results for the ERM (LRM). The isocontours correspond 
to probability of 5%, 10%. and 40%. The black points indicates the position in the (Wmaic, F°^^) plane of GRB 050904. The point with arrow means that the 
gap size should be considered as an upper limit, since the corresponding dark region could be affected by the presence of a DLA (Kawai et al. 2006). 




A direct comparison of the model with data can be carried 
on in terms of the probabiUty to find the largest gap in a given 
width range for burst afterglows at z = zgrb- When applied to 
the only known GRB at z > 6, i.e. GRB 050904 at zgkb = 6.29, 
a clear indication is obtained that reionization must have occurred 
well before z — 6. 

We find that the observed LGW in the GRB 050904 afterglow 
spectrum are consistent with xhi = 6.4±0.3 x 10~^. This result is 
in agreement with previous measurements by Totani et al. (2006), 
who find that a::Hi is bound to be xhi < 0.17 (< 0.6) at 68% (95%) 
C.L. 

Some gaps/peaks in absorption spectra could be due to 
DLAs/HIl regions intervening along the LOS towards the back- 
ground source (e.g. the DLAs/transverse proximity effect detec- 
tion by Totani et al. 2006/Gallerani et al. 2007). Such contam- 
inants could affect neutral hydrogen measurements, if the fore- 
ground sources are not properly removed. 

In this work we have applied the gap statistics to the GRB 
050904 observed flux. However, the same analysis can be done by 
using the transmitted flux Fobs/ F{v), i.e. normalizing the observed 
flux to the continuum. Following this approach, we provide useful 
plots (Fig. 2) which allow a straightforward comparison between 
our model results and future observations of 2 > 6 GRBs afterglow 
spectra. 

It is worth noting that the gap statistics with a variable flux 
threshold can be readily applied to other background sources, as 
QSOs. We plan to apply this analysis to the Fan et al. (2006) data, 
thus improving the results obtained by Gallerani et al. (2007). The 
advantage of varying the flux threshold is that this technique allows 
to study the response of the transmitted flux both to the high and 
low end tail of the IGM density field. In fact, high (low) Fth values 
probe regions corresponding to low (high) overdensities. 

In spite of the fact that QSO data are presently more abundant, 
GRBs are expected to be found at higher redshift, as they should re- 
sult from the death of (early) massive stars (e.g. Abel et al. 2002, 
Schneider et al. 2002). Larger samples of GRB afterglow spectra at 
z > 6, likely available in the near future, will allow a statistically 
significant analysis and possibly to reconstruct the cosmic reion- 
ization history. At redshifts approaching reionization epoch, gaps 
in afterglow absorption spectra are expected to become as large as 
the spectral region between the Lya and Ly/3 emission lines, thus 
obscuring completely the GRBs optical counterpart. 
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